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Abstract 

We discuss the production of photon pairs in hadronic collisions, from fixed target to LHC 
energies. The study which follows is based on a QCD calculation at full next-to-leading order 
accuracy, including single and double fragmentation contributions, and implemented in the form 
of a general purpose computer program of "partonic event generator" type. To illustrate the 
possibilities of this code, we present the comparison with observables measured by the WA70 
and DO collaborations, and some predictions for the irreducible background to the search of 
Higgs bosons at LHC in the channel h 77. We also discuss theoretical scale uncertainties for 
these predictions, and examine several infrared sensitive situations which deserve further study. 
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1 Introduction 



The production of pairs of direct photons|^ with large invariant mass is the so called irreducible 
background for the search of the Higgs boson in the two photon decay channel in the intermediate 
mass range 80 GeV < m/j < 140 GeV at the forthcoming LHC. This background is huge and 
requires to be understood and quantitatively evaluated. 

Beside this important motivation, this process deserves interest by its own. The production of 
such pairs of photons has been experimentally studied in a large domain of energies, from fixed 
targets |l], |2|, |3| to colliders @, ^ ^. A wide variety of observables has been measured, such as 
distributions of invariant mass, azimuthal angle and transverse momentum of the pairs of photons, 
inclusive transverse momentum distributions of each photon, which offer the opportunity to test 
our understanding of this process. 

The aim of this article is to present a study of diphoton hadroproduction based on a computer 
code of partonic event generator type. In this code, we account for all contributing processes 
consistently at next-to-leading order (NLO) accuracy, together with the so called box contribution 
99 ~^ 77- This code is flexible enough to accommodate various kinematic or calorimetric cuts. 
Especially, it allows to compute cross sections for both inclusive and isolated direct photon pairs, 
for any infrared and collinear safe isolation criterion which can be implemented at the partonic 
level. This article is organized according to the following outline. In section 2, we remind the 
basic theoretical ingredients, and present the method used to build the computer code developed 
for this study. Section 3 is dedicated to the phenomenology of photon pair production. We start 
with a comparison with fixed target and collider experiments. We then provide some predictions 
for LHC, together with a discussion of theoretical scale uncertainties. The theoretical discussion 
about the present day limitations of our code is continued in section 4. There we mention various 
infrared sensitive situations, which would deserve some more care, and for which the resummation 
of multiple soft gluon effects would be required, in order to improve the ability of our code to 
account for such observables. Section 5 gathers our conclusions and perspectives. 

2 Theoretical content and presentation of the method 

Let us first remind briefly the theoretical level of accuracy and limitations of works prior to the 
present one, in order to assess the improvements which we introduce. Then we present the method 
which we used to build our computer code DIPHOX. 

2.1 Theoretical content 

The theoretical understanding of this process relies on NLO calculations, initiated in Q. The 
leading order contribution to diphoton reactions is given by the Born level process qq — > 77 see for 
instance Diagram a. The computation of NLO contributions to it yields 0{as) corrections coming 
from the subprocesses qq — > 775, gq (or q) — > 775 (or q) and corresponding virtual corrections, see 

^The word "direct" means here that these photons do not result from the decay of tt", t), lu at large transverse 
momentum. Direct photons may be produced according to two possible mechanisms: either they take part directly 
to the hard subprocess, or they result from the fragmentation of partons themselves produced at high transverse 
momentum in the subprocess; see sect. 2. 
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for example Diagrams b and c. 




Diagram c 

Yet it also yields the leading order contribution of single fragmentation type (sometimes called 
"Bremsstrahlung contribution"), in which one of the photons comes from the collinear fragmen- 
tation of a hard parton produced in the short distance subprocess, see for example Diagram d. 
Prom a physical point of view such a photon is most probably accompanied by hadrons. Prom a 
technical point of view, a final state quark-photon collinear singularity appears in the calculation 
of the contribution from the subprocess gq 779- At higher orders, final state multiple collinear 
singularities appear in any subprocess where a high px parton (quark or gluon) undergoes a cascade 
of successive collinear splittings ending up with a quark-photon splitting. These singularities are 
factorized to all orders in as according to the factorization property, and absorbed into quark and 
gluon fragmentation functions to a photon D^/q ^j. g(z, Mj) defined in some arbitrary fragmenta- 
tion scheme, at some arbitrary fragmentation scale Mf. When the fragmentation scale Mf, chosen 
of the order of the hard scale of the subprocess, is large compared to any typical hadronic scale 
~ 1 GeV, these functions behave roughly as a/as{Mj). Then a power counting argument tells that 
these contributions are asymptotically of the same order in as as the Born term qq — > 77. What 
is more, given the high gluon luminosity at LHC, the gq (or q) initiated contribution involving 
one photon from fragmentation even dominates the inclusive production rate in the invariant mass 
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range 80 GeV < m^y < 140 GeV. A consistent treatment of diphoton production at NLO thus 
requires that 0{as) corrections to these contributions be calculated also, see for example Diagrams 
e and f. They have not been incorporated in 0, ^ and we compute them in the present work. 




+ ••• 



Diagram f 



The calculation of these corrections in their turn yields the leading order contribution of yet 
another mechanism, of double fragmentation type, see for example Diagram g. In the latter case, 
both photons result from the collinear fragmentation of a hard parton. In order to present a study 
of consistent NLO accuracy, NLO corrections to this double fragmentation contribution, see for 
example Diagrams h and i, have to be calculated accordingly. This is also done in the present 
article. 
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Diagram i 

We call "two direct" the contribution given by the Born term plus the fraction of the higher 
order corrections from which final state collinear singularities have been subtracted according to 
the MS factorization scheme. We call "one fragmentation" ("two fragmentation") the contribu- 
tion involving one single fragmentation function (two fragmentation functions) of a parton into a 
photon. Let us add one more comment about the splitting into these three mechanisms. One must 
keep in mind that this distinction is schematic and ambiguous. We remind that it comes techni- 
cally from the appearance of final state collinear singularities, which are factorized and absorbed 
into fragmentation functions at some arbitrary fragmentation scaleQ Mj . Each of the contributions 
associated with these three mechanisms thus depends on this arbitrary scale. This dependence on 
Mf cancels only in the sum of the three, so that this sum only is a physical observable. More 
precisely, a calculation of these contributions beyond leading order is required to obtain a (partial) 
cancellation of the dependence on Mf. Indeed this cancellation starts to occur between the higher 
order of the "two direct" contribution and the leading order of the "one fragmentation" term, and 
similarly between the "one-" and "two fragmentation" components respectively. This is actually 

^More generally, the definition of the fragmentation functions rely on the choice of a given factorization scheme, 
e.g. the MS scheme in this work. The fragmentation functions which we use are presented in fid]. 
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one of the first motivations of the present work. Thus, even though it may be suggestive to compare 
the respective sizes and shapes of the separate contributions for a given choice of scale, as will be 
done in p.2.1 , we emphasize that only their sum is meaningful. 




Diagram j 



Beyond this, the 0{a'^) so-called box contribution gg — > 77 through a quark loop is also in- 
cluded, see for example Diagram j. Strictly speaking it is a NNLO contribution from the point 
of view of power counting. However in the range of interest at LHC for the search of the Higgs 
boson, the gluon luminosity is so large compared with the quark and antiquark one, that it nearly 
compensates the extra powers of Og, so as to yield a contribution comparable with the Born term. 
For this reason, it has been included in previous works, and will be in the present one as well. We 
define the "direct" contribution as the sum "two direct" + box. 



Actually one should notice, firstly, that other NNLO gluon-gluon initiated processes, such as 
the collinear finite part of gg 9977 have been ignored^ although they could also be large. 
Secondly one should also even worry about the next correction to the box, because the latter may 
be quite sizeable. Such a possibility is suggested by the situation occurring to the first correction 
to the effective vertex gg h, computed in and shown to reach generically about 50 % of 
the one-loop result. Moreover, this box contribution is the leading order of a new mechanism, 
whose spurious (factorization and renormalization) scale dependences are monotonic, and only 
higher order corrections would partly cure this problem and provide a quantitative estimate. This 
tremendous effort has not been carried out yet, although progresses towards this goal have been 



achieved recently U2, QSl 14]. 



2.2 Presentation of the method 

In 0, a dedicated calculation was required for each observable. Since then more versatile ap- 
proaches have been developed, which combine analytical and Monte-Carlo integration techniques 
[^], |]l5|. They thus allow the computation of several observables within the same calculation, at 
NLO accuracy, together with the incorporation of selection/isolation cuts at the partonic level in 

■^The collinear divergent parts of these 2^4 processes have been already taken into account in the NLO corrections 
to the "one fragmentation" contribution and leading order "two fragmentation" components respectively. 
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order to match the various cuts used by the experimental collaborations as faithfully as possible. 
The studies of § and of Q rely on such an approach. Let us briefly describe the one which we use 
here. 



2.2.1 Phase space slicing and subtraction of long distance singularities 

Within the combined analytical and Monte-Carlo approach, two generic well known methods can 
be used to deal with infrared and collinear singularities which are met in the calculation of inclusive 
cross sections: the phase space slicing method [16] and the subtraction method [^]. The approach 
followed in the present work uses a modified version of the one presented in [15|, which combines 
these two techniques. 



For a generic reaction 1 + 2— >3 + 4 + 5 two particles of the final state, say 3 and 4, have a 
high pt and are well separated in phase space, while the last one, say 5, can be soft, or collinear 
to either of the four others. The phase space is sliced using two arbitrary, unphysical parameters 
PTm and R in the following way: 

- Part I 

The norm pT5 of transverse momentum of the particle 5 is required to be less than some 
arbitrary value pTm taken to be small compared to the other transverse momenta. This 
cylinder supplies the infrared, and initial state collinear singularities. It also yields a small 
fraction of the final state collinear singularities. 

- Part II a 

The transverse momentum vector of the particle 5 is required to have a norm larger than pTm, 
and to belong to a cone C3 about the direction of particle 3, defined by (ys — 1/3)^ + {(p^ — 
^s)^ ^ Rfh^ with Rth some small arbitrary number. C3 contains the final state collinear 
singularities appearing when 5 is collinear to 3. 

- Part II b 

The transverse momentum vector of the particle 5 is required to have a norm larger than 
PTm, and to belong to a cone C4 about the direction of particle 4, defined by (ys — 2/4)^ + 
((/)5 — (pi)'^ < R\^. C4 contains the final state collinear singularities appearing when 5 is 
collinear to 4. 

- Part II c 

The transverse momentum vector of the particle 5 is required to have a norm larger than 
PTm, and to belong to neither of the two cones C3, C4. This slice yields no divergence, and 
can thus be treated directly in 4 dimensions. 

Collinear and soft singularities appear when integration over the kinematic variables (trans- 
verse momentum, rapidity and azimuthal angles) of the particle 5 is performed on parts I, II a and 
II b. They are first regularized by dimensional continuation from 4 to d = 4 — 2e, e < 0. The 
d-dimensional integration over the particle 5 on these phase space slices yields these singularities 
as 1/e poles together with non singular terms as e ^ 0. After combination with the correspond- 
ing virtual contributions, the infrared singularities cancel, and the remaining collinear singularities 
which do not cancel are factorized and absorbed in parton distribution or fragmentation functions. 
The resulting quantities correspond to pseudo cross sections where the hard partons are unresolved 
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from the soft or collinear parton 5, which has been "integrated out" inclusively on the parts I, 
II a, II b. The word "pseudo" means that they are not genuine cross sections, as they are not 
positive in general. They are split into two kinds. We call pseudo cross section for some 2 — > 2 
process the sum of the lowest order term plus the fraction of the corresponding virtual corrections 
where the infrared and collinear singularities have been subtracted, and which have the kinematics 
of a genuine 2 — > 2 process. The contributions where the uncanceled collinear singularities are 
absorbed into parton distribution (on part I) or fragmentation (on parts II a and II b) functions 
involve an extra convolution over a variable of collinear splitting, as compared to the kinematics 
of a genuine 2^2 process: we call them pseudo cross sections for quasi 2 — > 2 processes. The 
detailed content of these terms is given in the Appendix^^. For an extended presentation of the 
details and corresponding explicit formulas, we refer to |[l5|| . 



As a matter of principle, observables do not depend on the unphysical parameters pTm and 
Rth- Yet, the pseudo cross sections on parts I, II a, II b and II c separately do. Let us briefly 
discuss the cancellation of the pTm and Rth dependences in observables computed according to this 
method. In the cylindrical part I, the finite terms produced are approximated in order to collect 
all the terms depending logarithmically on pTm, whereas terms proportional to powers of p^m are 
neglected. This differs from the subtraction method implemented in the cylinder in |T^, which 
kept the exact pTm dependence. On the other hand, in the conical parts II a and II b, the same 



subtraction method as in |15] is used, so that the exact Rth dependence is kept. This ensures the 
exact cancellation of the dependence on the unphysical parameter Rth between part II c and parts 
II a, II b whereas only an approximated cancellation of the unphysical parameter pTm dependence 
between parts II c, II a and II b and part I occurs. The parameter pTm must be chosen small 
enough with respect to pT3 and pn in order that the neglected terms can be safely dropped out. In 
practice, it has been verified that pTm values of the order of half a percent of the minimum of pTs 
and PT4: fulfill these requirements. A more detailed discussion on this issue is provided in Appendix 



The pseudo cross sections on parts I, II a, II b, as well as the transition matrix elements on 
the part II c, are then used to sample unweighted kinematic configurations, in the framework of a 
partonic event generator, described in 2.2.2 below. 



2.2.2 Partonic event generator 

For practical purposes, a partonic event generator has been built for diphoton production including 
all the mechanisms: the "direct", "one-" and "two fragmentation". Each mechanism is treated sep- 
arately. Firstly, the contribution of a given mechanism to the integrated cross section is calculated 



with the integration package BASES |18]. At this stage, some kinematic cuts (e.g. on the rapidity 
of the two photons, on their transverse momenta, etc.) may be already taken into account. Then, 
for the 2 — > 2 contributions and the quasi 2^2 contributions on the parts I, II a and II b, and 
the inelastic contributions on the part II c of the phase space, partonic events are generated with 
SPRING 1 18] with a weight ±1 depending on the sign of the integrand at this point of the phase 



spac^. All the events are subsequently stored into a NTUPLE fl^ . Finally these NTUPLES can 
be histogramed at will, incorporating any further cuts, such as those imposed by some isolation 

^This trick circumvents the fact that SPRING works only with positive integrands, while the pseudo cross sections 
are not positive. The generated events are thus unweighted up to a sign. 
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criterion as discussed in the next subsection. It is suitable to use values for Rth and pTm which are 
fairly small and disconnected from any physical parameter. The phase space generation is then as 
exclusive as possible. Moreover it allows to investigate the dependence of various observables with 
respect to the physical isolation parameters, as well as to investigate different types of isolation 
criteria, using an event sample conveniently generated once for all. In practice however one cannot 
use too small values in order to keep statistical fluctuations under control, unless the computer 
time and the sizes of the NTUPLES become intractably large. 



Let us state more clearly what we mean by partonic event generator. Since the events associated 
to the 2 — > 2 and quasi 2^2 contributions have a negative weight, this code, properly speaking, is 
not a genuine event generator on an event-by-event basis. By events, we mean final state partonic 
configurations. For a given event, the informations stored into the NTUPLE are the 4-momenta of 
the outgoing particles; their flavors: parton (i.e. quark or gluon) or photon; in the fragmentation 
cases, the longitudinal fragmentation variable(s) associated with the photon(s) from fragmentation; 
and, for practical purpose, a label which identifles the type of pseudo cross section (2 — > 2, quasi 
2 — > 2, inelastic) which produced the event stored. Notice also that in the fragmentation cases, all 
but the longitudinal information on the kinematics of the residue of the collinear fragmentation is 
lost. Hence this type of program does not provide a realistic, exclusive portrait of final states as 
given by genuine, full event generators like PYTHIA |2C] or HERWIG [21|. On the other hand, the 
latter are only of some improved leading logarithmic accuracy. Thus, our code is more precisely a 
general purpose computer program of Monte-Carlo type, whose virtue is the computation of various 
inclusive enough observables within the same calculation, at NLO accuracy. 



2.3 The implementation of isolation cuts 

Collider experiments at SppS, the Tevatron, and the forthcoming LHC do not measure inclusive 
photons. Indeed, the inclusive production rates of high prp vr", t], uj, or of pairs vr^Tr*^ or 77r'', etc, with 
large invariant mass, are orders of magnitudes larger than for direct photons. In order to reject 
the huge background of secondary photons produced in the decays of these mesons, the experi- 
mental event selection of direct photons (single photons, as well as diphotons) requires the use of 
isolation cuts. Such a requirement will be absolutely crucial at LHC for the search of Higgs bosons 
in the two photon channel and the mass range 90-140 GeV, since the expected background from 
tt", etc. is about eight orders of magnitudes larger than the signal before any isolation cut is applied. 

A widely used criterion to isolate photons is schematically the following^. A photon is said to 
be isolated if, inside a cone centered around the photon direction in the rapidity and azimuthal 
angle plane, the amount of hadronic transverse energy El^'^ deposited is smaller than some value 
ETmax fixed by the experiment: 

{y - y.,f + {(t) - <l,^f^ < \ 

Erp°- < Exmax J 

The topic of the isolation of photons based on the above cone criterion (|l|) has been rather 
extensively discussed in the theoretical literature, especially in the case of production of single 

^An alternative to the criterion (|l]) has been recently proposed in in which the veto on accompanying hadronic 
transverse energy is the more severe, the closer the corresponding hadron to the photon direction. It has been designed 
to make the "fragmentation" contribution vanish in an infrared safe way. 
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isolated photons in hadronic^ collisions ||2^, 24, 25, 26, |2^. Beside the rejection of the background 
of secondary photons, the isolation requirement also reduces the photons from fragmentation. The 
account of isolation effects on the "fragmentation" contribution was accurate to LO accuracy in 
1 23, A treatment to NLO accuracy has been subsequently given in [p^ j, following the subtrac- 
tion framework presented in [^5[. Isolation implies however that one is not dealing with inclusive 
quantities anymore. This raised questions concerning the validity of the factorization property in 
this case, and whether the fragmentation functions may depend on the isolation parameters, as 
assumed in [25|. This raised also issues regarding soft gluon divergences in isolated photons cross 
sections, as in |29]. These questions have been clarified in |27, |3^. The factorization property of 
collinear singularities still holds for cross sections based on the criterionj^ and the fragmentation 
functions involved there are the same as in the inclusive case, whereas the effects of isolation are 
consistently taken into account in the short distance part. Yet cross sections defined with this cri- 
terion may have infrared divergences - or, at least, instabilities, depending on the inclusiveness of 
the observable considered - located at some isolated critical points inside the physical spectrum of 
some observables calculated at fixed order, namely NLO, accuracy. This means that the vicinities 
of these critical points are sensitive to multiple soft gluon effects, which have to be properly taken 
into account in order to provide correct predictions. 



In the present calculation, as in ||2^, 27 1 for the case of single photon production, the transverse 
energy deposited in the cone may come from the residue of the fragmentation, from the parton 
5 (which never fragments into photons) or from both. During the projection of the NTUPLES 
onto any desired observable, the isolation criterion (|l|) about the two photons is applied to each 
stored partonic configuration. The effects of isolation are commented in 3.2.2. In addition, at the 



NLO accuracy at which our calculation is performed, potentially large logarithmic contributions 
of infrared origin may be induced by the extra isolation constraint on the phase space. The issue 
of infrared sensitivity induced by isolation will be discussed further in Let us mention that no 
summation of such logarithms is performed in our treatment. 



3 Phenomenology 

In this section, we adopt a LHC oriented presentation. We start with a brief comparison of our NLO 
calculations with WA70 and DO data for illustrative purposes. We then show some predictions for 
LHC in the invariant mass range 80 GeV< m^^ < 140 GeV corresponding to the Higgs boson search 
through /i — > 77. We discuss the ambiguities plaguing these predictions due to the arbitrariness in 
the choices of the renormalization scale fi, of the initial state factorization scale M (which enters 
in the parton distribution functions), and of the fragmentation scale Mj. 



3.1 Comparison with experimental data 



The related topic of isolated prompt photons produced in e e~ annihilation into hadrons has also been abundantly 
discussed. A variant of the definition (Q) suitable for e^e~ has been studied in j28|, and recently revisited in p9| , |3C| ]. 
An alternative criterion has been proposed in [Q, and applied to the measurement of isolated photons in LEP 
experiments. 

'^The fact that transverse energies are involved in in hadronic collisions is crucial in this respect. Factorization 
would be broken if energies were used instead. 
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3.1.1 Fixed target data 



A comparison between the diphoton differential cross section versus each photon's transverse mo- 
mentum measured by the WA70 collaboration ||l| and our NLO postdiction is shown on Fig. |l[ 
together with the respective magnitude of the various contributions. The NLO calculation has 
been made with the ABFOW parton distribution functions for the proton and the correspond- 
ing ones for the pion |33|^, for the scale choice^ M = Mf = fJ- = A(pt(7i) + Pr(72))! with A = 0.275. 
The "one fragmentation" contribution is one order of magnitude below the "two direct" contribu- 
tion. The "two fragmentation" contribution is even smaller and negligible here. The smallness of 
these contributions is the reason why previous works [0, ^ described this observable reasonably 
well too, despite the absence of higher order corrections to the fragmentation contributions there. 



Various correlations between the two photons: the distribution of the px imbalance variable 
z = — Pt(7i)-Pt(72)/Pt(7i) distribution of the azimuthal angle between the two photons {(l)yy), 
the distribution of PoMff^, and the distribution of transverse momentum of diphotons (qx), have 
been measured also by the WA70 collaboration [^. These distributions are infrared sensitive near 
the elastic boundary of the spectrum (e.g. — *■ or cp^^ — > vr) or near a critical point (e.g. 
z = 1) and, moreover, are quite sensitive to non perturbative effects appearing in the resummed 
part of calculations summing soft gluon effects. This sensitivity extends over a wide part of the 
spectrum covered by the measurements. Consequently we do not present any comparison of these 
data points with the approximation of fixed order accuracy of this work; nor will we discuss the 
scale ambiguities at fixed target energies. 



3.1.2 Tevatron collider data 

A preliminary study of diphotons events in the central region (|2/(7i,2)| < 1-0) has been recently 
performed by the DO collaboration Q . 

The experimental cuts in the DO data used for the comparisons are not corrected for electromagnetic 
calorimeter absolute energy scale. The electromagnetic energy scale correction is given by [^: 

S(measured) = a E(tiue) + 5 

*The choice of the parton distributions is mandated by the fact that the initial state of the reaction is vr* proton. 
Therefore a consistent set of parton densities inside the proton and the pion must be taken. Indeed, to extract 
the parton distribution functions in the pion, reactions such as tt^ p — > 7* (Drell-Yan) and tt* p ^ 7 X (direct 
photon) are used. Consequently some correlations between the proton and the pion partonic densities exist, and it is 
preferable to use consistent sets in the calculations. Only three groups provided such a work: ABFW [Q, MRS 
and GRV jssj. All these works are rather old and the partonic densities are rather similar in the WA70 x range. 

^One shall not attach importance to the somewhat unusual value A — 0.275 of the scale choice. Relatively low 
scales such as this one, or A = 0.25 equally well, turn out to match the data better than higher scale choices. Yet 
this particular value was not chosen as the one which matches the data the best, but for a minor though cumbersome 
computational reason. The WA70 collaboration requires the transverse momenta of the photons to be larger than 3 
GeV and 2.75 GeV respectively. However for computational convenience we first implemented a symmetric cut on 
the pt of each photon: pT > 2.75 GeV at the level of the Monte Carlo generation of photon pairs. In the ABFOW 
parametrizations, the factorization scale has to be larger than 2 GeV^. Given the above symmetric cut on both 
photons in the Monte Carlo generation, taking A = 1/4 does not ensure that is always above 2 GeV'^, while the 
choice A > .275 does. 

^'^The beam axis together with the direction of one of the two photons define a plane. The component of the 
transverse momentum of the other photon along the direction perpendicular to this plane is the pout of this photon. 
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where 

a = 0.9514 ± 0.0018l[J:[![]^^ 
5 = -0.158 ±0.015lJ]:^?GeV 

Thus, the experimental cuts at measured values of 14 (respectively 13) GeV correspond to cuts 
at roughly 14.90 (resp. 13.85) GeV in the theoretical calculation. Smearing effects accounting 
for electromagnetic calorimeter resolution have not been implemented, but given the experimental 
fractional energy resolution of the electromagnetic calorimeter p^ , they are expected to be of the 
level of a few percent only. 



The actual isolation cuts used experimentally (such as vetoes on charged tracks in some conical 
vicinity about each photon, etc.) are quite more complicated than the schematic criterion (||), and 
cannot be faithfully implemented at the partonic level. We instead simulated them in our NLO 
calculation by requiring that the accompanying transverse partonic energy be less than Exmax = 2 
GeV in a cone R = 0.4 about each photon. Varying Exmax from 1 to 3 GeV in the calculated 
cross-section, as a rough estimate of the effects of smearing due to hadronic calorimeter resolution 
and unfolding of underlying events contribution turns out to have a less than 4% effect. 

The MRST2 set of parton distributions functions^ is used^, with the scales arbitrarily 
chosen to be M = Mf = ^ = The prediction for the above scale choice is shown for 

the diphoton differential cross sections vs. the transverse momentum of each photon (Fig. ^), the 
diphoton mass (Fig. |3|), and for the transverse momentum of photon pairs (Fig. |^ and the az- 
imuthal angle between the photons (Fig. |5|). With the scale choice used, the "one fragmentation" 
contribution is roughly one tenth of the "direct" one whereas the "two fragmentation" yields a tiny 
contribution. To illustrate this, the different contributions: "direct", "one-" and "two fragmenta- 
tion" are shown separately on Fig. ^. The distributions of the transverse momentum of photon 
pairs and of the azimuthal angle between the photons are well known to be controlled by multiple 
soft gluon emission near the elastic boundary of the spectrum, qx ^ Q and vr respectively. 

Consequently, the accuracy of any fixed-order calculation, including the present one, is not suited 
to study such observables in these respective ranges. More on this issue will be commented in the 
next section. On the other hand a NLO calculation is expected to be predictive for the tails of 
these distributions away from the infrared sensitive region. 

The data are reasonably described, taking into account a correlated systematic error for events 
in which the px of both photons is above 20 GeV. This correlated systematic error due to the back- 
ground evaluation affects obviously the three highest pT points of the transverse energy spectrum, 
as well as the three highest points of the diphoton mass spectrum. 



^^The MRST sets 1,2,3 are associated with the value A-j^-? ~ 300 MeV for ti/ = 4 flavors. This corresponds to 
ctsinnz) ~ 0.1175 in the MS scheme. For more details, see [ p7[ 

^^The MRSTl set is presented by the authors of [B7| as the default set. However, in order to take into account 
mutually inconsistent data sets on single direct photon production at fixed targets, a kr smearing procedure is 
involved in the determination of this set. This procedure is strongly model dependent and questionable as long as 
no unambiguous way is found to lodge it in the QCD improved parton model. The set MRST2 does not involve this 
procedure, so we prefer to base any prediction and comparison on this set. 
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We do not present any analysis of the various scale dependences for Tevatron. Such a discus- 
sion is proposed for LHC in the next section. Yet let us mention that, at Tevatron, the energy 
scale is lower and the relevant values of x are somewhat higher than at LHC. Consequently, the 
renormalization scale dependence is slightly sharper, on the other hand the factorization scale de- 
pendence is somewhat flatter than at LHC. Nevertheless the situation at Tevatron is expected to 
be qualitatively similar to the one at LHC. 



3.2 Predictions for LHC 

We now discuss some results computed with the kinematic cuts from the CMS and ATLAS proposals 
0, namely pt{ii) > 40 GeV, pt{12) > 25 GeV, |y(7i,2)| < 2.5, with 80 GeV < m^^ < 140 GeV, 
and using the MRST2 set of parton distribution functions |37| and the fragmentation functions of 

0. 



3.2.1 Scale ambiguities 

We first consider the invariant mass distribution of diphotons, in absence of isolation cuts, cf. Fig. |6| 
in order to illustrate the strong dependence of the splitting into the three contributions, "direct", 
"one-" and "two fragmentation" , on the scale chosen, as we warned in In both choices of scales 
displayed the "one fragmentation" contribution dominates, but the hierarchy between "direct" and 
"two fragmentation" contributions is reversed from one choice to the other. With the choice of 
scales M = Mf = fj, = the "one fragmentation" is more than twice larger than the "di- 

rect" one, and the "two fragmentation" is the smallest. On the other hand, with the other choice 
M = Mf = fi = 2m^^, the "one fragmentation" contribution is three to five times larger than the 
"two fragmentation" component, and more than one order of magnitude above the "direct" one. 
On the other hand the total contribution seems rather stable. 



Yet the arbitrariness in the choices of the various scales still induces theoretical uncertainties 
in NLO calculations. In the following we actually do not perform a complete investigation of all 
three scale ambiguities independently with search for an optimal region of minimal sensitivity. At 
the present stage, we limit the study to an estimation of the pattern and magnitude of their effect 
on our results. We show how the scale ambiguities affect our prediction for the invariant mass 
distribution. We consider both the case without isolation (Fig. |^) and the isolated case with 
ETmax = 5 GeV inside R < 0.4 (Fig. ^). For the present purpose, the virtue of the actual values 
of the isolation parameters used here is to strongly suppress the fragmentation contributions hence 
the associated Mf dependence. We compare four different choices of scales: two choices along the 
first diagonal = M = Mf = and = M = Mf = Im^-y; and two anti diagonal choices, 

fi = m^^/2;M = Mf = 2m.y-y and n = 2m^^;M = Mf = m^^/2. We do not perform a separate 
study of fragmentation scale dependence. Yet the latter can be indirectly estimated by comparing 
the results of the isolated case, where the fragmentation components, thereby the corresponding 
fragmentation scale dependence, are strongly suppressed, with the situation in the non isolated 
case, where especially the "one fragmentation" contribution is quite large, and the "two fragmen- 
tation" not negligible, so that the issue of fragmentation scale dependence matters. 

When scales are varied between m^^/2 and 2m^^ along the first diagonal fi = M = Mf, the 
NLO results for the invariant mass distribution appear surprisingly stable, since they change by 
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about 5% only. Alternatively, anti-diagonal variations of /x and M = Mf in the same interval about 
the central value m^^ lead to a variation still rather large (up to 20 % cf. Fig. ^ and Fig. P). This 
is because variations with respect to ^ and M are separately monotonous but act in opposite ways. 
When n is increased, as{^'^) hence the NLO corrections decrease]^ On the other hand the relevant 
values of momentum fraction of incoming partons are small, ~ 0(10^'^ to 10~^) so that the gluon 
and sea quark distribution functions increase when M is increased. In the isolated case, this leads 
to a monotonous increase of the "direct" component, over a large band of the invariant mass range 
considered, as M is increased, cf. Fig. ^, which is induced in particular by the monotonous increase 
of the box contribution. Scale changes with respect to and M turn out to nearly cancel against 
each other along the first diagonal but add up in the other case. Actually, the stability along the 
first diagonal is accidental. 

In conclusion, the /i, M dependences are thus not completely under control yet at NLO in the 
kinematic range considered. On the opposite, the account for the NLO corrections to the fragmen- 
tation components provides some stability with respect to Mf variations about orthodox choices 
of the fragmentation scale. 



The issue of dependence of less inclusive observables, such as the tails of the qt or cj)^^ 
distributions are the same for the invariant mass distribution. This is because the tails of these 
distributions is purely given by the NLO corrections and dominated by the 0{as) corrections of the 
"two direct" component. On the other hand, the M dependence is a bit larger, so is the combined 
uncertainty on the theoretical results for these distributions, cf. Fig. O and Fig. 11. 



3.2.2 Effect of isolation 

We now consider the effect of isolation on the various contributions. As expected, isolation reduces 
the diphoton production rate, with respect to the inclusive case, cf. Fig. |l^. More precisely, severe 
isolation requirements like Exmax = 5 GeV inside a cone R = 0.4 suppress the "one fragmentation" 
component, which dominates the inclusive rate, by a factor 20 to 50, and kill the "two fragmenta- 
tion" contribution completely. 



However this net result hides a rather intricate mechanism, cf. Fig. 13 vs. Fig. g by which the 
"two direct" contribution turns out to be increased! Surprising as it may seem at first sight, this 
effect has the following origin. Higher order corrections to the "two direct" component involve in 
particular the two subprocesses qq 775 and gq — > 775 (where g is a quark or an antiquark) . The 
first one yields a positive contribution. On the other hand, the collinear safe part of the second 
one yields a contribution which is negative, and larger in absolute value than the previous one 
in the inclusive case, as was already seen in Isolation turns out to suppress more the higher 

In processes for which the lowest order is proportional to some power a", n > 1, an explicit jj, dependence appears 
in the next-to- leading order coefficient function, which partially compensates the (large) fj, dependence in Qs(/^^) 
weighting the lowest order. Unlike this, in the "two direct" component which dominates the cross section when a 
drastic isolation is required, the lowest order involves no as- This leads to a rather small fi dependence, since the latter 
starts only at NLO. On the other hand, the /i dependence occurs only through the monotonous decrease of the Qs(/i^) 
weighting the first higher order correction: there is no partial cancellation of fi dependence. Such cancellation would 
start only at 0(0^), i.e. at NNLO. The mechanism is more complicated in presence of fragmentation components, 
and the situation becomes mixed up between all components when the severity of isolation is reduced. 
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order corrections from the second mechanism than from the first one, so that the NLO isolated 
"two direct" contribution is larger than the inclusive one. Yet, the "fragmentation" contributions 
are suppressed more than the "two direct" one is increased, so that the sum of all contributions 
is indeed decreased, with respect to the inclusive case. Once again, one has to remember that the 
splitting into the three mechanisms depend, not only on the factorization scale, but more generally 
on the factorization scheme. This arbitrariness generates such counterintuitive offsprings; in a final 
state factorization scheme different than the MS scheme, the various components, especially the 
"two direct" one, may be separately affected by isolation cuts in a different way. This once more 
illustrates the danger of playing with these unphysical quantities separately. 

A more detailed analysis of the dependence of NLO estimations of various observables on the 
isolation cut parameters, especially on ETmax will be given in a forthcoming publication. We will 
also come back to this issue, regarding infrared sensitivity, in |4.3| . 



4 Infrared sensitive observables of photon pairs and soft gluon 
divergences. 

Being based on a fixed, finite order calculation, our computer code is not suited for the study of 
observables controlled by multiple soft gluon emission, and has to be improved in this direction. 
Among these infrared sensitive observables, one may distinguish the following examples, most of 
which would require an improved account of soft gluon effects. 



4.1 Infrared sensitivity near the elastic boundary 

4.1.1 The transverse momentum distribution da/dqT of photon pairs near = 

Both in the inclusive and isolated cases, this distribution is an infrared sensitive observable, con- 
trolled by the multiple emission of soft and collinear gluons. This well known phenomenon has been 



extensively studied for the corresponding observable in the Drell-Yan process [41|. A loss of bal- 
ance between the contribution of real emission, strongly suppressed near this exclusive phase space 
boundary, and the corresponding virtual contribution, results in large Sudakov-type logarithms of 
rr? jq^ (m being the invariant mass and and qr the transverse momentum of the photon pair - the 
heavy vector boson in the Drell-Yan case) at every order in perturbation. In order to make sensible 
predictions in this regime, these Sudakov-type logarithms have to be resummed to all orders. 



The treatment of the "two direct" and box contributions is similar to the well-known Drell-Yan 
process, and has been carried out recently by |42] at next-to- leading logarithmic accuracy in the 
framework tailored by Collins, Soper and Sterman |43]. On the other hand, the fragmentation 
contributions do not diverge order by order when — > 0. Indeed, in the "one fragmentation" case, 



parton^ + parton2 
parton3 



71 + parton3 

72 + ^ 



the NLO contribution to the hard subprocess (|2|) yields a double logarithm of the form 

~ as In^ ||Pr(7i) + PrCpartong)]! 



(2) 
(3) 



(4) 
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when ||pt(7i) + Pr(parton3)|| 0. However the extra convolution associated with the fragmenta- 
tion (^) involves an integration over the fragmentation variable pr(72)/pT(pa'rton3) which smears 
out this integrable singularity. The "two fragmentation" contribution involves two such convolu- 
tions, hence one more smearing. 



4.1.2 The distribution of photon-photon azimuthal angle da/d(j)^^ near (j)^^ = ir 

This distribution is another interesting infrared sensitive observable, measured by several exper- 
iments both at fixed target and collider energies |||, ^, ^, though less discussed in the literature 
from the theoretical side. The regime (p^yy — > vr includes back-to-back photons, a set of configu- 
rations which lie at the elastic boundary of the phase space. This case differs from the previous 
one for two reasons. Firstly, not only the "two direct" contribution diverges order by order when 
<j)^^ — > vr, but also both "one-" and "two fragmentation" contributions diverge as well, as can be 
inferred from Fig. |5[ Indeed, consider the example of the "one fragmentation case" , cf. equations 
§ and |3[ Selecting (p^^ — > vr emphasizes the configurations with (;/)(parton3) — 0(71) vr, so that 
all the emitted partons besides parton 3 have to be collinear to either of the incoming or outgoing 
particles, and/or soft, which yields double logarithms 

~ as In^ [vr - {(j){pavton^) - (5) 

associated with each of the hard partons 1,2,3 - plus single logarithms as well. For the observ- 
able da/d(j)'y'y near (p^^ = vr, the integral involved in the convolution of the hard subprocess with 
the fragmentation functions does not smear these logarithmic divergences, since the fragmentation 
variable pT(72)/PT(parton3) is decoupled from the azimuthal variable (j){psLrton^) which is equal to 
0(72)5 72 and parton 3 being collinear. A similar observation holds for the "two fragmentation" 
component. Moreover, in both fragmentation cases, soft gluons may couple to both initial and final 
state hard emitters. The resulting color structure of the emitters is more involved than in the "two 
direct" case, and especially more complicated in the "two fragmentation" case as shown in some 



recent works |44]. This would make any resummation quite intricate beyond leading logarithms. 



Let us notice that both fragmentation components make da/dcj)^^ diverge also when (p^^ — > 0. 
The increase of the fragmentation contributions in the lower (p^^ range is the trace of this divergence, 
cf. Fig. I 

4.2 An infrared divergence inside the physical region. 

In the case of photons isolated with the standard fixed cone size criterion of eqn. (|l]), a new problem 
appears in the qt distribution. This problem does not concern the region qt 0; still it has to do 



with infrared and collinear divergences. This can be seen on Fig. 14, which shows the observable 
da/dqT vs. qx for isolated photon pairs, computed at NLO accuracy. The computed qt distribution 
turns out to diverge when qt Exmax from below. Notice that the critical point Exmax is located 



inside the physical region. The phenomenon is similar to the one discovered in [25| in the production 
of isolated photons in e^e~ annihilation, and whose physical explanation has been given in 
following the general framework of |^] . It is a straightforward exercise to see that the lowest order 
"one fragmentation" contribution has a stepwise behaviour, as noticed in Indeed, at this order, 
the two photons are back-to-back. ^ad being the transverse hadronic energy deposited in the 
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cone about the photon from fragmentation, the conservation of transverse momentum imphes at 
this order that Exhad = Qt- The corresponding contribution to the differential cross section da/dqT 
thus takes the schematic form: 

r \ (1 /'"igm, LO) 

=f{qT)e{ET^ma.-qT) (6) 



According to the general analysis of |45|, the NLO correction to da/dqT has a double logarithmic 
divergence at the critical point qr = Exmax |3- The details of this infrared structure are very 
sensitive to the kinematic constraints and the observable considered. In the case at hand, at NLO, 
da/dqT gets a double logarithm below the critical point, which is produced by the convolution of the 
lowest order stepwise term above, with the probability distribution for emitting a soft and collinear 
gluon: 



\dqTj 



(1 fragm, HO) 

~ -/ (qt) {Exmax - qr) 



where C is a color factor, Cp or Nc according to whether the soft collinear gluon emitter is a quark 
(antiquark) or a gluon. More generally, at each order in Og, up to two powers of such logarithms 
will appear, making any fixed order calculation diverge at qT = Exmax, so that the spectrum com- 
puted by any fixed order calculation is unreliable in the vicinity of this critical value. An all order 
resummation has to be carried out if possible in order to restore any predictability. A correlated 
step appears also in the "two direct" contribution at NLO, in the bin about qx = Exmax- A 
detailed study of these infrared divergences will be presented in a future article. 



No such divergence appears in the qx distribution of photon pairs presented in The non 
appearance of the double logarithmic divergence there comes from the fact that the latter pops out 
only at NLO, while the authors of |Q compute the "one fragmentation" component at lowest order. 
Furthermore, the stepwise lowest order "one fragmentation" contribution to the qT distribution is 
replaced in by the result of the Monte Carlo simulation of this component using PYTHIA [20|. 
A quantitative comparison is thus difficult to perform^ 



In practice, the qr spectrum is sampled into bins of finite size, and the distribution represented on Fig. is 
averaged on each bin. Since the logarithmic singularity is integrable, no divergence is actually produced. However 
when the bin size is shrunk, the double logarithmic branch appears again. 

^^Such a comparison involves two issues. 
The first aspect concerns the infrared sensitivity below the critical point. When the scale of Qs in the Sudakov 
factor of the fragmenting quark is chosen to be the transverse momentum of the emitted gluon with respect to the 
emitter, the parton shower not only reproduces the fragmentation function of a parton into a photon to the collinear 
leading logarithmic approximation, but it also provides an effective resummation of soft gluons effects to infrared and 
collinear leading logarithmic accuracy. (This would not be true if, instead, the scale of as in the Sudakov factor were 
the virtuality of the emitter). This ensures that the distribution does not diverge from below at the critical point, 
but rather tends to a finite limit. 

The second issue concerns the shape of the tail of this contribution above the critical point. Indeed, energy-momentum 
conservation at each branching makes the parton shower generate also contributions in the region qt > Et max , which 
is forbidden at lowest order. These contributions would be classified in a beyond leading order calculation as higher 
order corrections. Unlike in a fixed order calculation however, they provide only a partial account of such corrections, 
but to arbitrary high order. The accuracy of these terms is thus uneasy to characterize, and a quantitative comparison 
between PYTHIA and any fixed order calculation is difficult to perform. 
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It can be noticed that the divergence at qx = Exmax is not visible on Fig. ^. This is because 
in this case, the critical point Exmax in the qx spectrum where the theoretical calculation diverges 
is too close to the other singular point qx = 0, given the binning used. The two singularities 
contribute with opposite signs in these bins and a numerical compensation occurs, resulting in no 
sizeable effect. Yet the problem is only camouflaged. A similar smearing appears also at LHC 
energies for a stringent isolation cut, cf. Fig. 

4.3 Reliability of NLO calculations with stringent isolation cuts 

Let us add one more comment concerning NLO partonic predictions with very stringent isolation 
cuts. In such calculations, the isolation cuts act on the products of the hard subprocess only. On 
the other hand, in an actual LHC event, a cut 2.5 GeV inside a cone R = 0.3 

or 0.4 will be nearly saturated by underlying events and pile up. 

This means that such an isolation cut actually allows almost no transverse energy deposition 
from the actual hadronic products of the hard process itself. This may be most suitable experi- 
mentally, and one may think about simulating such an effect safely in an NLO partonic calculation 
by using an effective transverse energy cut much more severe than the one experimentally used. 
However, requiring that no transverse energy be deposited in a cone of fixed size about a photon 
is not infrared safe, i.e. it would yield a divergent result order by order in perturbation theory. 
This implies that NLO partonic calculations implemented with finite but very stringent isolation 
cuts in a cone of fixed finite size would lead to unreliable results, plagued by infrared instabilities 
involving large logarithms of Ex max- What is more, these infrared nasties would not be located at 
some isolated point in the diphoton spectrum (like some elastic boundary or some critical point, 
as in the previous subsection), but instead they would extend over its totality, even for observables 
such as the invariant mass distribution. The issue of an all order summation of these logarithms 
of Ex max would have to be investigated in this case. 

5 Conclusions and perspectives 

We presented an analysis of photon pair production with high invariant mass in hadronic collisions, 
based on a perturbative QCD calculation of full NLO accuracy. The latter is implemented in the 
form of a Monte Carlo computer programme of partonic event generator type, DIPHOX. The post- 
dictions of this study are in reasonable agreement with both WA70 fixed target, and preliminary 
DO collider data, in the kinematical range where the NLO approximation is safe, namely away from 
the elastic boundary of phase space. Yet more will be learnt from the final analysis of the Tevatron 
data, and even more so after the Tevatron run II in the perspective of the LHC. It will then be 
worthwhile to perform a more complete phenomenological study. 

This notwithstanding, there remains room for improvements. A first improvement will be to 
take into account multiple soft gluon effects in order to calculate infrared sensitive observables 
correctly. Another improvement will concern a more accurate account of contributions beyond 
NLO, associated namely with the gluon-gluon initiated subprocess. Among those are the NNLO 
corrections, and even the two loop, so-called double box correction to gg — > 77, which may be 
quantitatively important at LHC for the background to Higgs search. 
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A better understanding of the effects of isolation, and their interplays with infrared problems 
is also required. This concerns the qt distribution iiGctr the criticctl point qj^ — E'j' rnax 

induced by 

isolation even when E^max is not small; this concerns also the status of partonic predictions when 
ETmax is chosen very small. Alternatively it would be interesting to explore the properties of differ- 
ent isolation criteria, such as, for example, the one invented recently by Prixione [^]. Concerning 
these last two items, approaches relying on beyond leading order partonic level calculations, and 
full event generators like PYTHIA or HERWIG will be complementary. 
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A Technical details on the two photon production 

In this appendix, we give some details on the method used to deal with infrared and soft divergences. 
For a complete presentation, we refer to [^]. The most complicated kinematics happens in the 
two fragmentation mechanism. Only the two fragmentation contribution will be treated in this 
appendix, the kinematics of the other cases can be simply deduced replacing the fragmentation 
function by a Dirac distribution: 



D^,^{x,Mj) = 5{l 



X 



At the hadronic level, the reaction Hi(Ki) + H2{K2) 7(^3) +7(i^4) is considered with: 

Ki = ^(1,0,1) 

K2 = ^(1,0,-1) 

K3 = Kt3 (coshys, ns, sinhya) 

K4 = Kt4 (cosh?/4, n4, sinhy4) 



where 

rig = n| = 1 

The cross section of the preceding reaction is the sum of the following parts. 

- The part I (cf. sect. 2.2.1) contains the infrared, the initial state, and a part of the final state 



collinear singularities. Once these divergences have been subtracted, i.e. cancelled against 
virtual divergences or absorbed into the bare parton distribution (for the initial state collinear 
singularities) or the bare fragmentation functions(for the final state collinear singularities), 
this part generates three types of finite terms. 

(i) The first type, of infrared origin, has the same kinematics as the lowest order (LO) terms 



and is given in A. 4 Pseudo cross section for the infrared and virtual parts. 
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(ii) The second type, of initial state collinear origin, has an extra integration over the center 
of mass energy of the hard scattering, as compared to LO kinematics. For this reason, 
it is called quasi 2 — > 2. It is given in A. 2 Pseudo cross sections for the initial state 
collinear parts. 

(iii) There is also a third type, of final state collinear origin, which involves also an extra 
integration as compared to LO kinematics. 

The parts II a and II b contain the rest of the final state collinear singularities. Once these 
divergences have been absorbed into the bare fragmentation functions, the remaining finite 
terms involve an extra integration over the relative momentum of the collinear partons, as 
compared to LO kinematics. These terms are combined with those of the so called third type 
(iii) above, cf. equations ( |A.10 ) and ( |A.ll ). The resulting contributions are called quasi 
2 ^ 2 as well. They are given in A. 3 Pseudo cross sections for the final state collinear parts. 



The part II c has no divergences. It is given in |A.l| Cross section for real emission. 



A.l Cross section for real emission 

The cross section is parametrized in the following way: 



O" 



Cij / dys / dy4 / dKrs / dKr^ 



dXZ fPTHmax 

/ dpT5 PT5 



PTr, 



^35— Cs 



X PT3 D^/k{x3,M 



D^/i{x4,Mf) 



PT4 

F,/H,{xi,M^) Fj/H^{X2,M^) 



X 



IMI 



Xl X2 

^ dX4 fPT5ma^ 

+ I — / dpT5PT5 

^4inin ^ PTm o /"* 



X 



PT4 D^/i{x4,Mj:) 



PT3 

2^ Z7 ('^_ i\/r2\ 



Fi/H^{xi,M') Fj/H,{X2,M' J 

X I 1 4 



Xl 



X2 



^'35 dy5 



'M5 dy5 



(A.l) 



where 
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2Ktz 

xsmin = ^ coshya (A.4) 
= ,A.5) 



s 

The transverse momenta pxs (resp. pta) are the transverse momenta of the fragmenting partons. 
They are related to the photon variables by pTs = Kts/x^ (resp. pn = Kxi/xi). The integration 
range for the pair of variables (p^^ (resp. ^45), 1/5 is the kinematically allowed range minus a cone 
in rapidity azimuthal angle C3 (resp. C4) along the P3 (resp. ) direction whose size is Rth- The 
overall factor Cij reads: 

and the Cj are given by: 

for quarks 



{N — 1) for gluons 



The matrix element squared taken from the first reference of and [^], has been split into 
two parts: 

|Mp = \M\l + \M\l 

The first part |-/Vf|| contains final state collinear singularities arising when / / P5 and the second 
part I Mil contains final state collinear singularities arising when p^ // p^. More precisely, the 
matrix element squared can be written as a weighted sum of eikonal factors Eat plus a term free 
of infrared or collinear singularities: 

4 4 

\M\l^Mn. = Yl Yl HabiP5)Eab + G{p5) (A.6) 
a=l b=a+l 

where 

Pa-Pb 
^ab — 



Pa-P5 Pb-Pb 

Using: 



1 1/11, , 

+ ----] (A-7) 



P3-P5P4-P5 Pl-Pb + P2-Pb \PZ-Pb Pi-Pb 

we get: 



^^An overall factor of the matrix element squared containing the average on spins and colors of the initial state 
and the coupling constant has been put into the coefficient dj 
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\M\l 



\M\ 



^Hi2{p5)Ei2 + i/i3(p5)^13 + H23{P5)E23 + ^34(^5)^34 + ^^(ps) 
^Hi2{p5)Ei2 + Hi4p5)Eu + H24{P5)E2A + ^34(^5)^34 + ^^(ps) 



with 



E. 



P3-P4 



1 



34 



-^34 



P1-P5 +P2-P5 P3-P5 
PI-P5 +P2-P5 P4-P5 



In order that the infrared divergences cancel, and the cohinear singularities factorize out, the 
coefficients Hab have to fulfill: 



Zl 



Huiil - zi)pi) + His({l - zi)pi) + HiJ{l - zi)pi 



^a^]{z2)\M\^P^,, 



Z2 



Hul^il- Z2)P2) + H23(^il - Z2) P2) +H2i({l- Z2)P2 



Z3 



Hn ( - — —P3] + H23 (- — —P3] +H34 (- — —p3 



Z3 



Z3 



Z3 



aS?(z4)|M|2f^,, 



Z4 



Hl4 { —P4] + H24 (- —P4] +H34 (- —P4 



Z4 



Z4 



Z4 



In particular, the cancellation of infrared divergences is insured by: 



^12(0)+/7l3(0) + ii'i4(0) 
^^12(0)+i/23(0) + i/24(0) 
^13(0)+i^23(0) + if34(0) 
//l4(0)+i/24(0) + i^34(0) 



\2B 



4?(1)|M|^,, 
4'^)(l)|M|?%. 



The functions af.\z) will be given in equation ( [A. 13 ) 



In equation ( [A.l ), the integration domain for the rapidities and the transverse momenta of the 
two photons is in general limited by experiments. The integration over px^ is constrained by: 



pl^<S{l-Xl){l-X2) 
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A. 2 Pseudo cross sections for the initial state collinear parts 

The finite part associated to the colhnear divergence Vx 1 1 given by: 



J J J J X^jYiin "^3 X^fjiin '^'^ 

X C| PT 6{pT3-pT4)D.,/k{x3,M]) D.,/i{x^,M]] 



27r 

x9 zi X? 



i'j^kl 



Cj/ 



(i-^i)+ 



+ In 



PXrn 

M2 



where the variables x\ (resp. x^) are defined by: 



Xi 



™0 

X2 



PT 



-S/3 



S 



and pt stands for pxz or pta- 

The finite part associated to the collinear divergence P2 // P5 is given by: 



^5//2 



dx^ 



X2. 



dxA 



^4n 



X4 



2tt 



-C§ PT d{pT3 - Pta) D^/k{x3,M]) D^/i{xA, M 



X 







^0 
X2 



ij'^kl 



£1 



<'j (^2) ^ 

+ In 



with 



1 - Z2) 



c. 



PTm 



M2 

4 52QC7^. 



(A. 



(A. 



The functio ns g.^ '^\z), pj^j\z) and will be defined at the end of this appendix cf. 

equations from (|A.18D to (|A.21| ), ( |A.13D and (1A.24|) . 



A. 3 Pseudo cross section for the final state collinear parts 

These parts contain the collinear singularities which have been absorbed into the bare fragmenta- 
tion functions. 
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The finite part associated to the collinear divergence // is given by: 

u u u u u X'trn.in 

Z3 



X O -/ \^^-^\ij^k'l 

Xi ^2 

whereas the finite part associated to the colhnear divergence P4 // P5 is given by: 

dX4 dZ4 



(A.IO) 



+ 2 



'^5//4 



J dyz j dy4 j dKrs j dKr^ j 



X4 



Z4„ 



Z4 



X Cf^ D^/,ixlM-f) D^,i{x4,Mj) 
X ~i 177 W^Hj^kl' 



+ 2 



•^1 -^2 
4'^~'^(^4) (pl4\ i4) 

In(l - Z4)\ (4) 



P^v^iz4)-dn'{z4) 



(A.ll) 



(1 - Z4) 



a];>{z4) + ln(i?2) e(z4™ - Z4) 



(1 - Z4) 



The functions a'^'f(z) and dij{z) will be also defined at the end of the appendix cf. equations 
from (lAll) to ( |A17|) and ( [Al25|) . The variables 

^3m ) ^4m j ^Jj^ , X2 and X2 appearing in equations 

( [A.lOp and ( |A.lip are given by: 



■23m 
Z4jYi 



PT3 



X 



I 

X2 



Xl 



X2 



PT3 + PTm 
PTA 

PT4 + PTm 
PT3 + PT5 

Vs 

PT3 + PT5 

Vs 

PT4 + PTb 

Vs 

PT4 + PT5 



(el/3 + gj/4) 

[e-y^ + e-y*] 



Vs 
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A. 4 Pseudo cross section for the infrared and virtual parts 

This pseudo cross section is given by: 



(Jir 



dys I dyA I dKT3 / dKr^ 



2 vr 



Cij / ~ / — PT 0{pT3 - PTi) 



X D^/k{x3,Mf) D^/i{x4,Mf) 



Fi,H,{xlM^) Fj,H,{xlM^) 



ln( ^\ {bkk + hu) + ) (6,, + h,,) 



\M\ 



2B 



Hi3{0) + Hu{0) + H23i0) + H2M + 2/734(0) 



_h^2(plr, 



+ 



4"' V S 
^34(0) 



^^13(0) + /7l4(0) + H23{^) + ^24(0) + 2/734(0) 



2/7i2(0) + i7i3(0) + i7l4(0) + ^23(0) + ^24(0) 



TT 



>i34(/)+^34(-y'') +F{s,t,u) 



(A.12) 



with pt = Pt3 = PT4- The terms ba are defined in equations ( A.22| ) and ( A.23| ). In the equation 
( A.12| ), y* = {y3 — yi)/2 and the function A{x) is given by: 



A{x) = TT ln(2) ln(4 cosh^(y'')) + 2y* sinh(2y*) 



ln(sin ( 



+ 4 / dcp 

Iq cosh(2 x) + cos(2 ( 



ln(sin</>) arctan 



cosh(2 x) + cos(2 (/)) 
sin 6 



1 — cos ( 



The function F is the finite part of the virtual term and the variables s, i and u are the Mandelstam 
variables of the 2 — > 2 processes: 

S = ipUp'2f 
i=(pl-pl? 

u={pl-pl? 

where the 4- vectors p^ are the infrared limits of the 4- vectors pi. 



A. 5 Altarelli-Parisi Kernels 

We will give in this appendix the expressions of the functions aij(z) and bij. These functions are 
defined by: 
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{1-z 

{d-i), ^ 

where P^^-''' (resp. Plf*) are the AltareUi-Parisi Kernels in four (resp. d) dimensions. So the 
expressions for the functions a^^^{z), a^^~^\z) and hij are given by: 



4J(z) = 2 7V + (A.14) 



a^fjiz) = Cf{1 + z^) (A.15) 



"qq 

^2 



= cA^-^^-^]{l-z) (A.16) 



aW(z) = Tp{z^ + {l-zf){l-z) (A.17) 

where is the number of colors, Cp = (A^^ — 1)/(2A^) and Tp = 1/2. The extra part needed to 
get the functions a in d dimensions {o!f^{z) = o!f^ {z) — ea^f- *\z)) is given by: 

af-^\z) = (A.18) 

ag-^)(z) = Cpil-zf (A.19) 

ag-^)(z) = Cpzil-z) (A.20) 

ag-^)(z) = 2TFzil-zf (A.21) 

The coefficients bij read: 

. ffi^^ (A.2) 

bqq = ^Cf (A.23) 

The function fij{x) and cZij(^;) define the factorisation scheme for respectively initial state and 
final state collinear singularities. In the MS scheme, we have: 

fijiz) = (A.24) 

dij{z) = (A.25) 



B Cancellation of the pTm and Rth dependences 

In this appendix, we give further details on the cancellation of the pTm and Rth dependences in 
observables calculated according to the method used in this article. 
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In the conical parts II a and II b, the d-dimensional integration over particle 5 in Ci, i = 3,4, 
reads schematically: 

rPTSmax 
■J PTm 

i4>i5 dyr, sm — — — (B.l) 

Q cosh (yj - 2/5) - cos 

The term generating the final state collinear pole {p^//Pi) has been explicitly written, and the 
remaining quantity F(pT5,(j)i5,y5) is a regular function. In the parts II a and II b, the same 
subtraction method as in [l^ is used, and the following contribution is added and subtracted: 

J PTm JC, [Vi - Vb) + (Pi5 

In the cylindrical part I, the finite terms produced by the integration over particle 5 are approx- 
imated: all the terms depending logarithmically on pTm are kept, whereas terms proportional to 
powers of pTm are neglected. Notice that this differs from the subtraction method implemented in 



the cylinder in ||15[, which kept the exact pTm dependence. 



In summary, the present method is an admixture of the phase space slicing and subtraction 
methods, at variance with what has been done in It ensures the exact cancellation of the 

unphysical parameter Rfh dependence between part II c and parts II a, II b whereas only an ap- 
proximated cancellation of the unphysical parameter pTm dependence between parts II c, II a and 
II b and part I occurs. 

We checked carefully that the dependences on the unphysical parameters drop out. This point is 
illustrated by the pTm dependence (at fixed Rth = 0.1) and the Rth dependence (at fixed pTm = 0.1 
GeV), of the higher order (HO) part of integrated cross section (the lowest order (LO) part being 
independent of these parameters) 



a"^ = I dm. 



rrimax ^^HO 



15| and 16. We display separately the qq and qg initiated contributions to the "di- 



shown on Figs. 

rect" on Fig. [l5|, and the "one-" and "two fragmentation" mechanisms, on Fig. [l^. To be definite, 
the integration bounds are taken to be rrimin = 80 GeV, mmax = 1500 GeV, the cuts pT3, Pta > 25 
GeV, |y3,4| < 2.5 are applied, and the MRST2 set of parton distribution functions with the scale 
choice M = n = Mj = are used; let us emphasize however that the pattern obtained does 

not depend on these details. 

The quantity a^'^ does not depend on Rth and, in principle, it becomes independent of pTm at 
small enough pTm- To show these features more clearly, the observable displayed is the ratio Rm^^ 
defined as follows: 

Rm,, = ^ / dm^7 -7-— (B-3) 

^ Jrn^in " '"'77 

The integrated cross section is normalized to be asymptotically 1 in order to show the size of 
the relative error bars. However taking the denominator A equal to the calculated a^'-' for the 
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smallest value of pTm may be numerically unsuitable. Indeed, when pTm becomes smaller and 
smaller, numerical cancellations between larger and larger contributions occur and the error bars 
coming from the Monte Carlo integration become larger and larger. These numerical fluctuations 
affect the behavior in the limit pTm — 0. In order to bypass these technical complications, A is 
taken to be the averaged value of those of the integrated cross sections da^'-' /dm^^ which are 
consistent with each other within the error bars. For instance, for the pTm dependence of the 
"direct" contribution, the average is taken over the values corresponding to the three smallest PTm 
because the fourth one is not consistent with the others in the error bars. In addition, in the 
case of the direct contribution, the two partonic reactions qq and qg have been split because, for 
the above choices of scales, the two integrated contributions are large and of opposite signs. As 
expected, Rm^-^ does not depend on Rth and approaches 1 as PTm 0. Let us notice that one 
can wonder whether large relative fluctuations do not appear again when the two contributions of 
the "direct" are added. Indeed, the relative fluctuations of the HO terms are larger for the sum 
than for each parts, but these HO terms are small compared to the LO part 0(1%) (7 ) 

hence the "physical" cross section (LO+HO) is sufficiently stable. When the parameter pTm is 
chosen small enough with respect to prz and pta, the neglected terms power behaved in pTm can 
be safely dropped out. In practice, we observe that prm values of the order of half a percent of 
the minimum pTz and Pta^ i-c PTm < 0.1 GcV, fulfill these requirements. Before embarking in a 
long phenomenological study, the user of the DIPHOX code is advised to check whether the value 
of the parameter pTm to be used is small enough to neglect safely the power corrections of pTm- 
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Figure 1: Diphoton differential cross section da/dpT vs. pT^ the transverse energy of each photon, 
in 7r~-proton colhsions at ^/S = 22.9 GeV. Data points from the WA70 collaboration The 
solid line is the full contribution with scales M = n = Mj = 0.275 {prili) +Pt(72))- 
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Figure 2: Diphoton differential cross section da/dpT vspT, the transverse energy of each photon, at 
Tevatron, y/S = 1.8 TeV. PreUminary data points (statistical errors and systematics in quadrature) 
from the DO cohaboration |^] are compared to the theoretical predictions: the full NLO prediction 
is shown as the solid line. The ratio data/(full NLO theory) is shown below. 
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Figure 3: Diphoton differential cross section da/dra^-y vs. m^-y, the mass of the photon pair, at 
Tevatron, ^/S = 1.8 TeV. Prehminary data points (statistical errors and systematics in quadrature) 
from the DO collaboration |^] are compared to the theoretical predictions: the full NLO prediction 
is shown as the solid line. 
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Figure 4: Diphoton differential cross section da/dqT vs. qt^ the transverse momentum of the 
photon pair, at Tevatron, ^/S = 1.8 TeV. Prehminary data points (statistical errors and systematics 
in quadrature) from the DO collaboration Q are compared to the theoretical predictions: the full 
NLO prediction is shown as the solid line 
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Figure 5: Diphoton differential cross section da / d<j)'y'y vs. (p^-y, the azimuthal angle between the two 
photons, at Tevatron, = 1.8 TeV. Preliminary data points (statistical errors and systematics 
in quadrature) from the DO collaboration Q are compared to the theoretical predictions: the full 
NLO prediction is shown as the solid line while open squares (open circles) represent the single 
(double) fragmentation contribution. 
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Figure 6: Splitting of the diphoton differential cross section da/drri'y^ at LHC, VS = 14 TeV 
without isolation, into the "direct" , "one fragmentation" and "two fragmentation" components, 
shown for two different choices of scales. The following kinematic cuts are applied: prisfi) > 40 
GeV, pt(72) > 25 GeV, |y(7i,2)| < 2.5. 
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Figure 7: Diphoton differential cross section da/dm^'y vs. m^'y, the invariant mass of photon pairs, 
at LHC, s/S = 14 TeV without isolation. The following kinematic cuts are applied: prili) > 40 
GeV, pt{12) > 25 GeV, |y(7i,2)| < 2.5. The scale dependence is shown on the bottom plot. 
M = Mf is understood. 
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Figure 8: Diphoton differential cross section da/dm^^ vs. m.^^ at LHC, \/5 = 14 TeV, with 
isolation criterion Exmax = 5 GeV in i? = 0.4. Same kinematic cuts as in fig. 0. The scale 
dependence is shown on the bottom plot. M = Mf is understood. 
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Figure 9: M dependence of the "direct +box" contribution to da/dm^^ in several m^^ bins at 
LHC, \fS = 14 TeV, with isolation criterion Exmax = 5 GeV in i? = 0.4. Same kinematic cuts as 
in fig. 1^ /i is chosen to be while M is varied between and 2m^^. 
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Figure 10: Diphoton differential cross section da/dgx at LHC, a/S" = 14 TeV, with isolation 
criterion Exmax = 5 GeV in i? = 0.4. The following kinematic cuts are applied: prili) > 40 GeV, 
^^(72) > 25 GeV, 1^(71,2)1 < 2.5, and 80 GeV < m^j < 140 GeV. The scale dependence is shown 
on the bottom plot. M = Mf is understood. 
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Figure 11: Diphoton differential cross section da/d<j)^'y vs. (^^.y, the azimuthal angle between the 
two photons, at LHC, ^/S = 14 TeV, with isolation criterion Exmax = 15 GeV in i? = 0.4. Same 
kinematic cuts as in fig. |10|. The scale dependence is shown on the bottom plot. M = Mf is 
understood. 
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Figure 12: Diphoton differential cross section da/dm^^ vs. m.y.y at LHC, \^ =14 TeV, without 
and with isolation criterion ETmax = 5 GeV in i? = 0.4. Same kinematic cuts as in fig. |. The 
scale choice is M = Mj = /U = 
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Figure 13: Splitting of the diphoton differential cross section cZa/dm-y-y at LHC, VS = 14 TeV 
with isolation criterion Exmax = 5 GeV in i? = 0.4, into the "direct" , "one fragmentation" and 
"two fragmentation" components, shown for the scale choice ji = M = Mf = The following 

kinematic cuts are applied: pt{ii) > 40 GeV, pt{12) > 25 GeV, |y(7i,2)| < 2.5. 
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Figure 14: Diphoton differential cross section da/dqT at LHC, = 14 TeV, with isolation 
criterion Exmax = 15 GeV in i? = 0.4. Same kinematic cuts as in fig. ^. 
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Figure 15: Dependence of the ratio Rm-y^ (see equation ( |B.3D ) over the phase space shcing param- 
eters Rth and pTm for the "direct" contribution. 
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Figure 16: Dependence of the ratio Rm^^ (see equation ( |B.3D ) over the phase space shcing param- 
eters Rth and pTm for the "one-" and "two fragmentation" contributions. 



45 



